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1. Introduction 

In a previous paper [9] we have studied the expansion of an analytic function at two 
finite points in the complex plane. The domain of convergence is a Cassini oval around 
the two points. The main motivation for that paper was to obtain the coefficients of 
asymptotic expansions of certain integrals. In the present paper we give a few examples 
in which the expansion of an integral at two saddle points yields a convergent expansion 
that has an asymptotic property for large values of a parameter. 

In the well-known methods for deriving asymptotic expansions of integrals a basic 
step is transforming the integral into a standard form, and the transformation usually 
gives a new integral in which the integrand contains implicitly defined functions that 
are difficult to handle. In the method of this paper we avoid a transformation and, in 
addition, we derive convergent expansions. 

We start with a simple example in which only one saddle point occurs, and in which 
a function is expanded at that saddle point. This gives an expansion for the Charlier 
polynomials. 
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In two other examples (Laguerre and Jacobi polynomials) we take into account two 
saddle points, and again two convergent expansions can be constructed with the desired 
property. The approximants belong to the same class of polynomials as the original ones, 
but they are of a simpler type (Hermite and Chebyshev, respectively). The asymptotic 
property follows from recursion relations for functions appearing in the expansions. The 
convergence follows from the fact that an integral along a finite contour is expanded 
inside a domain of uniform convergence. 

In the examples given in this paper the contour integrals are based on Cauchy-type 
integrals obtained from generating functions. When the contour is finite a proof of the 
convergence is usually rather easy. For more general finite contours and more general 
integrals we expect that the method can be applied as well. For example, we can apply 
the method to the Gauss hyper geometric function and the incomplete gamma function, 
with different integral representations. 

Also, the methods of this paper can be generalized by considering Taylor expansions 
at more than two points. In [10] we give details on the theory of multi-point Taylor 
expansions, and in a future paper we will give details on applications to integrals with, 
for example, three saddle points. 

We show a few graphs that indicate the nature of the approximations, and in a final 
section we mention a few examples in which other functions are considered. 

2. A simplified version of the saddle point method 

Throughout this paper we are concerned with finding asymptotic expansions of integrals 
of the form 



where f(w) and g(w) are analytic in a domain Q of the complex plane that contains the 
origin; T is a circle with center at the origin and contained in f2; n is a large positive 
integer. We assume, as it usually happens to be the case, that the asymptotic behavior 
of the integral F(n) for large n is determined by contributions from the saddle points 
of <p(w) = g(w) -\nw [[15], Chap. 2, §4]. 

The standard saddle point method consists of 

(i) deforming the contour of integration V into a new path that crosses one or some 
of the saddle points of <p(w); 

(ii) a suitable change of the variable of integration; 

(hi) application of Watson's lemma or Laplace's method. 

Instead of applying the standard saddle point method, we will proceed in a simpler 
way: just substitute a power series expansion at one or more saddle points of the function 
f(w) in (2.1). If there is just one saddle point wq, then that power series is its Taylor 




(2.1) 
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expansion at wq: 

/H = E Z ^(-"-o) fc , (2.2) 
k=o fc - 

which is uniformly convergent for w in a disk D r {wo) = {w G Q, \w — wo\ < r} with 
center at wq and radius r = \rd we ^^\w — wq\. If there are two saddle points w\ and 
W2, then that power series is its two-point Taylor series at w\ and W2 [9]: 

oo 

f(w) = K(w - Wl ) + a' n (w - w 2 )\ (w - Wl ) n (w - w 2 ) n , (2.3) 

n=0 

where 

_ f(w 2 ) , _ fjwx) 
a = , a = (2.4) 

W2 — W\ W\ — W2 

and, for n = 1, 2, 3, 

1 ^ (n + fc - 1)! (-l)" +1 n/("- fc )(w 2 ) + (-l) fc fc/("~ fc )( Wl ) 
° n ~ n! ^ fc!(n-fc)! {w x - w 2 ) n+k+1 ' l ' j 

, = 1 (n + fc - 1)! (-l)" +1 n/("- fc )(wi) + (-l) fc fc/("- fc )( W2 ) 
° n n!^ fe!(n-fc)! (w 2 - ^i)"+ fc + 1 ' l ' j 

The expansion (2.3) is uniformly convergent for w in a Cassini oval 

O r (wi,w 2 ) = {w£Q, \w — wi\\w — W2\ < r} 

with foci at u>i and W2 and "radius" r = lni w€ ^ n {\w — Wi\\w — w 2 |}; see [9]. 

If we substitute now (2.2) or (2.3) in (2.1) and interchange summation and integra- 
tion we obtain an expansion of F{n). This is proved in the following two propositions. 
Proposition 2.1. Let the right-hand side of (2.2) converge uniformly to f{w) for 
w G D r (w ) with \ w \ < r, then 

F(«) = V£^ f(w-w ) k e^^. (2.7) 

Proof. If \wo\ < r, then G D r (wo). Then we can choose a small enough circle T in (2.1) 
such that T G D r (wo). Therefore, expansion (2.2) is uniformly convergent for w G T. 
Introducing (2.2) in (2.1) and interchanging summation and integration we obtain (2.7). 
□ 

Proposition 2.2. Let the right-hand side of (2.3) converge uniformly to f{w) for 
w G O r (wi, W2) with l^iti^l < f, then 



00 „ , 
F(n)=Y / a k / (w-wi) k+ \w-w 2 ) k e nifi{w) — + 

00 „ 

yx / (w-wi) k (w - w 2 ) 



k -° (2.8) 

k + l c nip(w) aw 
W 

k=0 
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Proof. The proof is similar to that of Proposition 2.1. □ 
In the remaining part of the paper we apply Proposition 2.1 or Proposition 2.2 
to three specific examples of integrals F(n) representing Charlier polynomials C^(nx), 
Laguerre polynomials L"(nx) and Jacobi polynomials P^ ,l3 \x). In this way, we obtain 
expansions of these polynomials for large values of n. In each example separately we 
prove that the corresponding expansions (2.7) or (2.8) are convergent in a certain region 
of the variable x and that in fact, they have an asymptotic nature for large n, uniformly 
with respect to x in certain domains of the region of convergence. 



n=0 



3. Asymptotic expansions of Charlier polynomials in terms of Gamma func- 
tions 

The Charlier polynomials are defined by the generating function 

e— '(1 + «,)* = f;C^(x)^, (3.1) 
and have the explicit expression 

n 

^(-) = E(y(IJ fe! (-«r- fc - (3.2) 

The Charlier polynomials are orthogonal with respect to a discrete distribution on the 
positive real line. For an overview of properties, see [7]. Recent papers on asymptotics 
are [1], [2], [5], and [6]. In this section we give a simple convergent expansion of C^(xn), 
that has an asymptotic property for large n, uniformly for complex a in compact sets 
and for complex x bounded away from 1. A uniform expansion that holds for x in a 
compact neighbourhood of x = 1 is given in [1], where J— Bessel functions are used in 
the approximations. Also, uniform expansions that hold for — oo < x < oo are given in 
[2]- 

Theorem 3.1. For x ^ 1, a G (D and n G N, the Charlier polynomials have the 
expansion 

OO , N J, 

C a n (xn) = e^ 1 "*) -^T^(x, n), (3.3) 

k=0 

where 

^ / ^ T(nx + 1) $ (x,n) 

$o(s,n) = rf v <S>i(x,n) = - y 3.4 

T{nx + 1 — n) (1 — x)(n(x — 1) + 1) 

and, for k = 0, 1, 2, 

Y(nx + 1) 1 

$ fc (x,n) = — — -r 2 F 1 (-k,-n,nx - n + 1;1 - x). (3.5) 

T(nx — n + 1) (1 — x) K 
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The sequence {&k(x, n), k = 0, 1, 2, ...} satisfies the recurrence 

(3.6) 



$fc(x,n) = 



n(x — 1) + k 



x(l — k) — k ^ x(l — k) , . . 

■* fc _i(x,n) + t -^-$ fc _ 2 (x,n) 



x-1 ■ (x _ 1)2 

where k = 2, 3, . . . . and is an asymptotic sequence for large n. For fixed k we have 

$ fc (x,n) = 0(n-L( fc + 1 )/ 2 J( n |x|) n ) (3.7) 

when n — > oo, where \a\ is the integer part of the real number a. The asymptotic 
property holds uniformly with respect to complex x, \x — 1| > 5 > 0. 

Observe that in the notation of Pochhammer's symbol (a) n , defined by 

(a) = l, (a)„ = V ; = a(a + l)---(a + n-l), n = 0,l,2,..., (3.8) 
r(o) 

we have 

$ (x) = (-lr ^"""^ = (-l) n (-nx)„ = (nx)(ns - 1) • ■ ■ (nx - (n - 1)). (3.9) 
1 (—nx) 

Proof. From (3.1) we derive the integral representation 

n\ f ^ dw 



C a n {x) 



j e- aw {l + wY 



2ni J T w n+1 ' 

where T is a circle with center at the origin and radius < 1. We write this in the form 

C%(nx) = — f e- aw e n ^ x ' w) — , (3.10) 
2iri J r w 

where <p(x, w) = x log(l+w;) — logu> . The only saddle point of tp(x, w) is wo = (x — 1) _1 . 
The function e~ aw is an entire function of w. Hence, the expansion 

e ~ aW = £ ( w ~ w o) k (3.H) 

fc=o K - 

is locally uniformly convergent for w € <Z*and x ^ 1. Therefore, after substituting this 
expansion in (3.10) and using Proposition 2.1 we obtain (3.3) with 

= ^!j / r ^-^) fe (l + ^r^T- (3-12) 
To obtain the recurrence (3.6) we write 

27Ti X — 1 J r 



(n-1)! 1 ,jl j. 9 e n^(x, W ) 

' " : 1 w) dir. 

ow 
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0.000000090 
0.100006223 
0.200157621 
0.301812498 
0.410358953 



0.534449998 
0.680932968 
0.855641877 
1.068772397 
1.347867376 



Table 3.2. The zeros of C^(nx) for n = 10 and a = 1. 



Integrating by parts and performing a few straightforward manipulations we obtain 
(3.6). Equalities (3.4)-(3.5) follow after simple calculations. The asymptotic behavior 
in (3.7) for large n follows from (3.4) and (3.6). From (3.4) we see that $o = 0(( n \ x \)n) 
(see formula (19) below) and that <I>i = 0((n\x\) n /n). Therefore, (3.7) is true for 
k = 0, 1. From here, the proof follows by induction over k. If (3.7) holds up to k, then 
= 0(n-W 2 \n\x\) n ) and $ fc = 0(n-[( fc+1 )/ 2 l (n\x\) n ). Using this in (3.6) with k 
replaced by k + 1 we have that &k+i = 0(n~^ k+2 ^ 2 \n\x\) n ). □ 

Property (3.7) holds uniformly for \x — 1| > 5 > 0, also for complex x. Also, detailed 
information on the asymptotic behavior can easily be obtained from (3.12). □ 

From (3.9) it follows that &o(x) has n zeros at x = m/n, m = 0,1, . . . n — 1 ($i(a?) 
has the same zeros, except for x = (n — l)/n). See also Table 3.2, where we give 
numerical values of the zeros of C^(nx) for n = 10 and a = 1. From the graphs in 
Figure 3.3 we also see that the early zeros are approximated quite well. 




FIGURE 3.3. Numerical experiment on the approximation given in Theorem 3.1 for large n 
and x £ [0, 1). Continuous lines represent the Charlier polynomial C^inx) for (a) n = 20 and 
(b) n = 50. Dashed lines represent the first order approximation given 
Both graphics are cut for extreme values of the polynomials. 



Remark 3.4. When in expansion (3.11) the expansion point wo is not equal to the 
saddle point (x — 1) , we are not able to prove the asymptotic nature of expansion 
(3.3). This follows from the integration by parts procedure mentioned in the proof of 
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Theorem 3.1. On the other hand, we can show (3.7) directly from the definition (3.12) 
with a change of variable like in the standard saddle point method. 

3.1. Details on the convergence 

It is of interest to verify the speed of convergence of the expansion in (3.3). We 
consider (3.12), with k = nn, and determine the saddle point of (w — wq) k (1 + w) x w , 
where k is large. We consider x fixed, and for x < 1 we verify that a positive saddle 
point w + occurs with w + ~ 1/[(1 — x)k]. There is a negative saddle point, whch is not 
relevant. We have 

[( W+ - W0 m +w+ r w - + r- {1 _ x)k } ■ 

Multiplying this with a k n\/k\, see (3.3) and (3.12), and using in Stirling's approximation 
of the factorials only the dominant parts, that is k\ ~ (k/e) k , we see that the main 
information on a k n\&k(k,n)/k\ is given by 

e k k n a k (l -x) n 
(1 — x) k k k ' 

where k is large compared with n and x, x < 1. We see that the ratio of successive 
terms is about a /[(I — x)k\. 

For other values of x, also complex, a similar analysis can be given, with some care 
in choosing the saddle points and defining the branches in the complex plane. 



n 


C^(0.25n) 


C a n (x,0) 


C%(x,l) 


CS(x,2) 


C™ (x, 3) 


C»(x,4) 


C° (x, 5) 


10 


-1.03630 


-0.97736 


-1.02335 


-1.04747 


-1.00438 


-1.03633 


-1.0363 


30 


4.35872 


4.03823 


4.28867 


4.35077 


4.35762 


4.35858 


4.35870 


50 


-4.86727 


-4.65813 


-4.82829 


-4.86464 


-4.86701 


-4.86725 


-4.86726 


90 


-2.94851 


-2.87926 


-2.93699 


-2.94808 


-2.94848 


-2.94851 


-2.94851 



Table 3.5. Numerical experiment on the convergence rate of expansion (3.3) for x = .25 
and a = 1. Here, C^(x,iV) = e"^ 1 "^ Ylk=o ^~k\ ^k{x,n), represents the truncated 
series in (3.3). All the rows are multiplied by an appropriate constant in order to keep 
the numbers small. 



4. Asymptotic expansions of Laguerre polynomials in terms of Hermite poly- 
nomials 

The Laguerre polynomials can be defined by the generating function 

oo 

(l_ t )-«-i e -W(i-*) = J^L°(x)t n , a,xe(D, \t\<l. (4.1) 

n=0 
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and have the representation 

n / \ k 

«w-E(-«*f:!:)ir- < 4 - 2 > 



k=0 



For deriving the asymptotic expansion we use the Cauchy integral that follows from 
(4.1): 

4°)(x) = — / e W(--i)(i _ w )- a -iJ^ (4.3) 

where T is a circle around the origin with radius < 1. The many- valued functions 
(1 — w)^ appearing here and in the theorem assume the principal branch that is equal 
to 1 at w = 

The asymptotics for large n, fixed a, is considered in [4]. For real x two uniform 
expansions are given, one involving the J— Bessel function for x in an interval that 
contains the origin, and one in terms of the Airy function for x in an interval containing 
the transition near the largest zero of L°(x). In this section we give an asymptotic 

1 /2 

expansion of L"(nx) in terms of L n (nx), which in fact is an Hermite polynomial. We 
consider x > 1 and for these values the expansion is convergent and is in particular 
of interest because this is the interval that contains the large zeros and the transition 
point at x = 4. 

When the parameter a of the Laguerre polynomial is large the asymptotic behavior 
can be described in terms of Hermite polynomials (see [12]). For example, we have the 
limit 

(—~\ \ n 9 — n / 2 

lim a~ n/2 L*{x^ +a) = ± '— HJx/\f2). (4.4) 

In [8] we have extended this limit by giving an asymptotic representation for large a 
and n fixed in terms of Hermite polynomials. For more details on large a asymptotics 
we refer to [12]. The approach in this section is quite different because we take a fixed 
and n large. 

From (4.3) we obtain 

n<p(x,w) rfw 



Li a) {nx) = -L / /H^ (4.5) 

Z7TI J T (1 — W) 6 I Z W 



where 



<p(x,w) = -logw, f(w) = (l-w) 1/2 - a . (4.6) 



The function ip(x, w) has two conjugate saddle points: 

x . i 



w 



± 



1 "2 ± 2 C ' £=VM4^). (4.7) 

The square root defining £ is positive for < x < 4; for x > 4 we define £ = iy / x(x — 4), 
again with positive square root. In the expansion of the Laguerre polynomials we allow 
that the saddle points coalesce. 
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4.1. Construction of the expansion 

The function f(w) of (4.5) is analytic in Q, = (D\ [1, oo) and we can expand f(w) in 
in a two-point Taylor expansion at the two saddle points w ± , using a slightly different 
form of (2.3), 

oo 

f(w) = Y J [A k + B k w] (w - w+) k (w -w~) k . (4.8) 

k=0 

After substituting expansion (4.8) in (4.5) and interchanging summation and integration 
we obtain 



L£\xn) = Y,[Ak*k(x,n) + B k * k (x,n)], (4.9) 

k=0 

1 r p nip(x,w) A w 

,,(,,„) _ _ jj w _ w+)Hw _ (4. 10 ) 

^ k (x,n) = — {w-w + ) k {w-w-) k - vrndw. (4.11) 

2m Jy (1 — w) ,i i' L 



(—l) n 

$o(x,n) = 4 1/2) (n*) = n[22n+1 r^ H 2n+1 (V^) , 



where 



and 



We have 



( _ ir -i ( 4 - 12 ) 
*o(a;,n) = i„Y 2) (nx) = ————=H 2n -i {Vnx) , 

and, for k = 1, 2, 3, 

* fc (s,n) •) x*-iL^;f N, * fc (x,n) - £ ( ■) ^"^^iM- 

J'=0 ^ ' j=0 ^ ' 

(4.13) 

The sequences {<3?fc(x, n)} and {^(x, n)}, fc = 0, 1, 2, . . ., satisfy the recurrences 

$fc = n _ 2 fc + 3 /2 {ai*fc_i + a 2 $ k -2 + 6i* fc -i + 62^-2} , (4.14) 

Gl = (& - l)(x 2 - 2x - 2) - i, a 2 = (k - l)x(2 - x), 

h = (k - 1)(2 - 3x) + — b 2 = (k - l)x(Ax - x 2 - 2), 

= o, 1 , ; /0 {co$fc + ci$ fc _i + c 2 $fc- 2 + + d 2 * fc - 2 , } , (4.15) 

n — 2k + 1/2 
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c = (2 - 3k)x + 2{k - 1) + 

x — 1 

ci = (1 - £:)x 3 + 4(fc - l)x 2 + kx + 2(1 — fc) + — ^— , c 2 = -6 2 , 

d 1 = (4ife - 3)x 2 + 2(4 - 5fc)x + 2(fc - 1) + 

d 2 = (k- l)x(x 3 - 6x 2 + 9x - 2). 
To verify the recursions (4.14) and (4.15) we write 

$ fc (x,n) = / (w - w + ) k ~ 1 (w - w~) h ~ 1 Vl - w dw, (4.16) 

2m n J c ow 

I r Q e n-ip(x,w) 

\& k (x,n) = / (w - w + ) k ~ 1 (w - w~) k ~ 1 wVl - w dw, (4.17) 

2m n Jc Ow 

1 f e nip(x,w) ^ u 

* fc _i(x,n) -w + $ fe _i(x,n) = — / (w - w+) k (w - w'f- 1 - r^ — 

2m Jc (1 — w) 6 ' z w 



(4.18) 



and 



* fc _i(x,n) -w;-$ fe _i(x,n) = — / (to - w+f" 1 ^ - uT) 7; — • 

2m J c (1 — w) 6 ' 2, w 



(4.19) 



Integrating by parts in (4.16) and (4.17), using (4.18) and (4.19) and after straightfor- 
ward manipulations we obtain (4.14) and (4.15). Formulas (4.12)and (4.13) follow from 
(4.10) and (4.11) after simple calculations. 

Theorem 4.1. Expansion (4-9) is convergent, uniformly for a G C in compact sets, 
and x > 1 + S > 1. Moreover, {<3?fc(x,n)} and |\l/fc(x,n)} are asymptotic sequences for 
large n: 

* k (x,n) =0 (n-L( fc +D/ 2 J) [\* (x,n)\ + \V (x,n)\] , 

; ; (4.20) 

* fc (x,n) =0 („-L(*+D/2J j [|<D (x,n)| + |^ (x,n)|] , 
as n — > oo and A; = 0, 1, 2, .... 

Proof. We apply Proposition 2.2. Expansion (4.8) is uniformly convergent for w 
inside the Cassini oval with foci w + and w~ and "radius" r = — — w~\, 

where w\ = 1 is the singular point of /. Using (4.7) it follows that r = x. The points w 
are inside the Cassini oval if they satisfy \w — w + \\w — w~\ < r = x. Because w + w~ = 1, 
the origin w = is inside the oval only if x > 1. Hence, the contour T of (4.5) can 
be taken completely inside the oval only if x > 1 (see also Figure 4.2). This proves 
the convergence of (4.9) for x > 1. The asymptotic behavior in (4.20) follows from 
(4.12) and the recursions (4.14) and (4.15). More detailed asymptotic information can 
be obtained from the integrals in (4.10) and (4.11). □ 




x = 7.5, x= 8 (lemniscate) and x = 8.5 

FIGURE 4.2. Cassini ovals for the expansion (4.8) for several values of x. For w inside the 
ovals, the expansion is convergent. For < X < 1 the origin is outside the oval; for x = 4 it 
is a circle, for x = 8 a lemniscate. For x > 8 the oval splits up into two parts. All ovals go 
through the point w = 1, a singular point of f(w). 

A 

A 




(a) (b) 
FIGURE 4.3. Numerical experiments on the approximation of Theorem 4.2 for large n and 
x £ (0,4]. Contin uous lines represent the Laguerre polynomial Ln (nx) for (a) n = 10 and 
(b) n = 20. Dashed lines represent the first order approximation given by Aq<&q(x, n) + 
Bq^q(x, n). Both graphics are cut for extreme values of the polynomials. 
Remark 4.4. The expansion in (4.9) has a meaning for all complex x, and has for all 
fixed x an asymptotic meaning. The expansion is uniformly convergent for |ac| > 1 + 8 > 
1. 

4.2. Details on the coefficients 

The expressions (2.5) and (2.6) for the coefficients in (2.3) can be used in the present 
case also. We first write 

/H = S K( w - w ~) + a 'k( w - w+ )] K w -w-)(w- w + )] k , (4.21) 

k=0 
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and compare this with (4.8). By comparing coefficients of equal powers, it follows that 
A k and B k can be expressed in terms of a k and a' k . We have for k = 0, 1, 2, . . . 

A k = -a k w + - a' k w~, B k = a k + a' k . 

We have 

Ak = -2» (w+ofc) , 5 fc = 25tf(a fc ), a = i(l - ^-) 1/2 " Q r \ (4.22) 
and, for = 1, 2, 3, 



; (fc + j-l)!(«-l/2) fe _,- f (-l)ij (-l) fc A; 



^ fc!j!(fc-j)!(»O fc+J ' +1 l(l-U;+) a + fc -J- 1 /2 ( 1 _ w -)a+fc-j-l/2 

(4.23) 

The coefficients can also be computed from the recursion relations 
x{k + l)A k+1 -x{k + l)B k+1 = (a-l/2 + 2k)A k -(xk + l)B k , 

(4.24) 

x(k + l)A k+1 - x(x - 3)(fc + l)5 fc+ i = (a + 1/2 + 2fc)B fe , 

where fe = 0,1,2,... . Let, for 1 < x < 4, x = 4sin 2 (6>/2). Then w ± = e ±ie , and 

A = -2 P+1 woP{6/2) cos[(# - tt)/?/2 -5], £ = 2 /3+1 waP(0/2) cos[(# - tt)/?/2]. 

where /3 = 1/2 — a. This gives real expressions for the first coefficients to start the 
recursion relations in (4.24). For x > 4 we can obtain expressions in terms of hyperbolic 
functions by writing x = 4 cosh 2 (0/2), which gives = —e ±e and 

A = -2^ +1 coslrV/2) cosh[(fl(/3/2 - 1)], B = 2^ +1 cosh' 3 (0/2) cosh(6>/?/2). 



4.3. An alternative form of the expansion 

By using in (4.5) the substitution 

f(w) = a + [3 Q w + (w- w~)(w - w + )g (w), (4.25) 
where a and /3 follow from substituting w = , we obtain by integrating by parts 

L ( C\nx) = a^ {x 1 n) + (3^ {x 1 n) + —— / h{w)- 7^ — , (4.26) 

2mn J v (1 — w) 6 > z w 

where 

AW = w(l - wf' 2 -^ [{y/T^g^w)] . (4.27) 
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Continuing this procedure we obtain the expansion in negative powers of n: 



oo oo 



0k 



k ' 



(4.28) 



fc=0 fc=0 



where a k and follow from 



f k (w ) = a k + (3 k w , f k {w + ) = a k + k w 



(4.29) 



where with fo(w) = f(w) and for k = 0, 1, 2, . . . 



/fe(ty) = a fc + (w - to )(u> - u> + )si fc (u>) 



A+iH = w(l - w) 3/2 — [(VI- wfffcH] • 



(4.30) 



The expansion in (4.28) also follows from re-arranging expansion (4.9) by using the 
recursion relations for <& k (x,n) and ty k (x,n) in (4.14) and (4.15). 



n 


(3.5n) 


(3.5,0) 


(3.5,1) 


lL 1} (3.5,2) 


(3.5, 3) 


1^(3.5,4) 


Ll 1} (3.5,5) 


10 


0.340506 


0.343249 


0.341724 


0.340495 


0.340449 


0.340490 


0.340504 


30 


-8.94039 


-8.86531 


-9.03530 


-8.95798 


-8.94213 


-8.94045 


-8.94038 


50 


-5.05678 


-5.05941 


-5.06764 


-5.05801 


-5.05689 


-5.05680 


-5.05678 


90 


6.56556 


6.56328 


6.57572 


6.56601 


6.56547 


6.56553 


6.56556 



Table 4.5. Numerical experiment on the convergence rate of expansion (4.9) for x = 
3.5 and a = 1. Here, L^(x, JV) = J2k=o^ k ^ k ( x ' n ) ^k^k(x, n)], represents the 
truncated series in (4.9). All the rows are multiplied by an appropriate constant in 
order to keep the numbers small. 



5. Asymptotic expansions of Jacobi polynomials in terms of Chebyshev poly- 
nomials 

The large n asymptotics for the Jacobi polynomials is discussed in [[3], Vol. II, §10.14], 
in particular for x G (—1,1). For x bounded away from the points ±1 elementary 
functions (sine and cosine functions) can be used for describing the asymptotics. For x 
close to ±1 Bessel functions can be used (Hilb-type formulas). 

In this section we develop a convergent expansion that is valid for x G (— 1, 1) and 
the terms of the expansion constitute asymptotic scales for large n. It is possible to 
extend the results to complex values of x, but this will not be considered here. The first 
approximants are Chebyshev polynomials, which in fact are elementary functions, and 
the other terms can be obtained from recursions that show the asymptotic property. 
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5.1. Construction of the expansion 

Starting point is the integral representation that follows from [[3], Vol. II, p. 172]: 

p(*,f}) M _ 1 (-!)" f (l-w-x) a (l + w + xf nvix , w) dw 
n {> 2m 2™ J r {l-x)"(l + x)f> w' { ' 

where we consider x £ (—1, 1); the function ip{x,w) is defined by 

<p(x, w) = log(l + w + x) + log(l — w — x) — log w 

and T is a simple closed contour, in the positive sense, around w = 0. The points 
w = — x±l are outside the contour, and (1 — w — x) a /(l — x) a and (1 + w + x) 13 / '(1 + x) 13 
are to be taken as unity when w = 0. 

The function tp(x, w) has two conjugate saddle points: 



= ±wo, wo = i\J 1 — x 2 . (5-2) 
We expand the integral by using the function 

f(w) = (1 - x)- Q -5 (1 + x)- p -*{l - w - x) a+ ^{l + w + xf + ^ , 
which is analytic in 

(D\ {(-oo, -1 - x] \J[1 - x, oo)}. 
We expand, using a slightly different form of (2.3), 

oo 

/H = l A k + B k w] (w 2 - wlf (5.3) 

fc=0 

where the coefficients A k and can be expressed in terms of the derivatives of f(w) 
at w = ±wo; see the next subsection for more details. 

After substituting expansion (5.3) in (5.1) and interchanging summation and inte- 
gration we obtain 

oo 

P^ P) ix) = [Ak*k(x, n) + B k V k (x, n)] , (5.4) 

k=0 

where 

ZJx n) - (-^^^ f (™ 2 + l-* 2 ) k c n M dw 

$fc(X ' nj " 2m 2" J r W(x,w) 6 w' (5 ' 5j 

y k ( x , n ) = L^l i [ (w 2 + l-x 2 ) k W(x,w)e n ^ x ^ — , (5.6) 

2m 2 n Vl - x 2 Jr w 
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and 

W(x,w) = ^ (I - w - x){\ + w + x). (5.7) 

We have 

d> (x,n) = Pt 1/2 '- 1/2) (x), * (x,n) = - x 2 )P^ 1,2 \x). (5.8) 

These Jacobi polynomials are Chebyshev polynomials: 

pt i/2 >- i/2) (*) = 2 -^$^T n (x), p(w/*) (a0 = 2 ;7 ( ( n 2 ;+; )! ^ ( x). (5 . 9) 

In terms of elementary functions: 

T n (cos#) = cosntf, U n (co S 9) = sm ( n + 1 ) g . ( 5 . 10 ) 

Furthermore 

1 - x 2 

$i(x,n) = [4(n + l)<f> (x,n) - Ax(n - l)V (x, n) + (2n - l)* (x,n- 1)] , 



(5.11) 

}> 

(5.12) 



_3(1 _ x 2\ 

*i(rc, n) = + 'i)( w + 2) ^ n + 1 ~ 2x n ]*o(x, n) + x(2n - l)* (x, n - 1)} , 
and, for k = 2, 3, 4, 

i=o v J 7 

»*<*.») - -(l-^tO) <l^|>^^^'«(x,. (5.14) 

j=o v J 7 

The sequences {<3?fc(:c, n)} and (\l/ fc (x,n)} satisfy the recursion relations 

$fc = n + 2 k _ l i a ^k-i + a 2 $k-2 + 6i*fc-i + &2*fc-2] , (5.15) 
= — zt [co^fc + ci$ fc _i + c 2 $ fc _ 2 + di* fc _i + d 2 ^fc- 2 ] , (5.16) 

where 

oi = (1 - x 2 )(6k - 5), a 2 = -4(fc - 1)(1 - x 2 ) 2 , 

= - x (4fc - 3), 6 2 = 4x(ife - 1)(1 - x 2 ), 
co = -x(4k-l), c 1 =x(l-x 2 )(8k-5), c 2 = -4x(l - x 2 ) 2 (£; - 1), 
di = 3(2k - 1)(1 - x 2 ), d 2 = -4(1 - x 2 ) 2 (fc - 1). 
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For the relations between Jacobi and Chebyshev polynomials we refer to [[13] p. 152- 
153]. The expressions in (5.11) and (5.12) follow from contiguous relations of the Jacobi 
polynomials (see[[13] p. 166]). 

To verify the recursions in (5.15) and (5.16), we write 

— 1 (— l) n r Qprnp{x,w) 
$ k ( x , n) = v - ' / (w 2 + 1 - x 2 ) fc_1 W(x, w) dw, (5.17) 

* fc (x, n) = 1 - ; / w(w 2 + 1 - x 2 )^ 1 W(x, w) — -= dw, (5.18) 

Integrating by parts in (5.17) and (5.18) and after straightforward manipulations we 
obtain (5.15) and (5.16). The asymptotic behavior pointed out above follows from the 
definition of the Jacobi polynomials and the recurrences (5.15) and (5.16). 

Theorem 5.1. Expansion (5.4) is convergent for x £ (—1,1). Moreover, {^(x, n)} 
and {^k(x,n)} are asymptotic sequences for large n: 

$ fc (x,n) = O (n-L( fc +D/ 2 J) , V k (x,n) = O (n~ ^ k+1 ^ ) , 

as n — ► oo and k = 0, 1, 2, .... 

Proof. Expansion (5.3) is uniformly convergent for w inside the Cassini oval of focus 
w^ and "radius" 

r = min{|l + x + w + \\l + x + w~\, |1 — x + w + ||l — x + w~\} = 2(1 — |x|). 

Hence, (5.3) is convergent for w G (D such that |u> 2 +l— x 2 | < 2(1— \x\). From Proposition 
2.2 it follows that expansion (5.4) is convergent for 1 — x 2 = \w + w~\ < 2(1 — |x|), that 
is, V x G (—1, 1). The asymptotic property follows from the recursion relations in (5.15) 
and (5.16). □ 

5.2. Details on the coefficients 

The expressions for (2.5) and (2.6) for the coefficients in (2.3) can be used in the 
present case also. We first write 

oo 

fH = J2 " + a 'k( w + w o)] ( w2 - w l) k i ( 5 - 19 ) 

k=o 

and compare this with (5.3). By comparing coefficients of equal powers, it follows that 
Ak and Bk can be expressed in terms of and a' k . We have for k = 0, 1, 2, . . . 

A k = (a' k - a k )w , B k = a k + a' k . 
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After straightforward manipulations we find 
29 [a k Vl 



A, 



where 



ar 



(1 _ x )a+l/2(! + x )/3+l/2 



and Bh = 



25fta fc 



(l_ x )a + l/2( 1+2; )/3+l/2' 

(5.20) 



a 



fl i \a+l/2 n . \8+l/2 

,l-o; + u>o) (l + x-u>o) 
2^n 



and, for m = 1, 2, 3, 



m ^ — ' 



(m + k- 1)1 



ml kl(m - k)l(2w ) m + k+1 



m — k 

£ 

3=0 



m — k 
j 



k{ _ ir+J{ _ a _ 1/2) ^_ (3 _ 1/2) 



m — k—j 



'1 — X — Wo 



.j-a-1/2 



(1 + X + Wq 



.m-k-j-0-1/2 



m (-l)*+J(-a- 1/2), (-0-1/2) 



m — k—j 



/i i \j~ a ~ 1/2 t-\ i \m—k—j— 3— 1/2 [ 

(1 — X + U>o) (1+X — Wo) J 



(5.21) 



(5.22) 




(a) (b) 
FIGURE 5.2. Numerical experiments about the approximation of Theorem 5.1 for large n and 
— 1 < x < 1. Continuous lines represent the Jacobi polynomial Pn (x) for (a) n = 10 and 
(b) n = 20. Dashed lines represent the first order approximation given by Aq^q + Bq^q. 



n 


p(3/2,l/2) (0) 


Pn(0,0) 


P»(0,1) 


Pn(0,2) 


P„(0,3) 


P»(0,4) 


Pn(0,5) 


10 


-0.336376 


-0.348029 


-0.348029 


-0.337153 


-0.336376 


-0.336340 


-0.336360 


30 


-0.201847 


-0.204304 


-0.204304 


-0.201909 


-0.201839 


-0.201845 


-0.201847 


50 


-0.157618 


-0.158781 


-0.158781 


-0.157636 


-0.157615 


-0.157617 


-0.157618 


90 


-0.118124 


-0.118612 


-0.118612 


-0.118128 


-0.118123 


-0.118124 


-0.118124 



Table 5.3. Numerical experiment on the convergence rate of expansion (5.4) for x = 0, 
a = 3/2 and (3 = 1/2. Here, P n (x,N) = Y^k=o\^ k ^ k ( x > n) + Bk^k(x, n)], represents 
the truncated series in (5.4). 
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6. A few other examples of convergent asymptotic expansions 

Consider the integral for the modified Bessel function (for properties of the special 
functions in this section we refer to [3] or [13]) 

{2z)"e* f\_ 2ztr 



U*) = rv ~T\ / e- M [t{l-t)\"-*dt. (6.1) 

The usual method for obtaining the asymptotic expansion for large z consists of 
substituting the expansion (1 — t) v_ 2 = ^ {"~ k *)(— t) h and by interchanging the order 
of summation and integration. The resulting integrals are not evaluated over [0, 1] but 
over [0,oo). This latter step gives the divergent asymptotic expansions. If we integrate 
over [0, 1] we obtain an expansion in terms of incomplete gamma functions: 

lv{z) - v^r(, + i)^ ( 1} { k ) (2z)* • (6 ' 2) 

For fixed values of z the incomplete gamma functions have the asymptotic behaviour 
j(v + k + 1/2, 2z) = y ' 1 [1 + Oik- 1 )] , k - oo, (6.3) 

and we see that the terms behave like 0(k~ u ~ 3 ^ 2 ). So, convergence is guaranteed if 
%tv > — \ . A further examination of the terms of the expansion shows that for large 
z it is better to use the asymptotic property of the expansion (the ratio of successive 
terms is of order 0(l/z)). The incomplete gamma functions can be computed by using 
a backward recusion scheme. 

Another example is the K— Bessel function given by 



K v {z)= v ; / e-^[t{l + t)Y-*dt. (6.4) 

I>+2) JO 

Again expanding (t+l)^ -1 / 2 = ^ (" ~ 5 ) (t) k gives the standard expansion. A convergent 
expansion can be obtained by using 

(t + ir-V 2 = £ Cfe (_A_)\ Ck = {-l)*(*- V y (6.5) 
This gives a convergent expansion in terms of confluent hypergeometric functions 



V k=0 



(6.6) 
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We have 

u(0,\-v,2z)=l, U(l, 1 --v,2z)=(2zy + ?e 2z r(- 1 --v,2z), (6.7) 

again an incomplete gamma function. Other U— functions can be obtained by recursion. 
For large k the [/—function behaves like [11] 

k\ U{k,\- v, 2z) = O (k a e - 2V *™^ , 

where a is some constant. It follows that the convergence is better than in the previous 
example. 

As a final example we consider an expansion of the Kummer function. Tricomi 
[14] has derived several convergent expansions of the i-Fi— function in terms of Bessel 
functions that are useful for evaluating the function when the parameters are large. For 
example, we have [14] 

oo n/2 

1 F 1 {a,c-z)=e^T(c){K Z ) {1 -^/ 2 Y J M^\c) (j£f J c _i+„ (2y^) , (6.8) 

n=0 

where k = \c — a and the A n (n, A) are coefficients in the generating function 

oo 

e 2KZ (1 - z) K ~ x (1 + z)- K ~ x = An(n, A) z n . (6.9) 

n=0 

The series in (6.8) is convergent in the entire z— plane. Moreover, it can be used for the 
evaluation of \Fi(a, c; z) for large k, because the series has an asymptotic property. For 
further details on these expansions we refer to [14]. 
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